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In the quasi-twodimensional (Q2D) electron gas of an InAs channel between an AlSb sub- 
strate and superconducting Niobium layers the proximity effect induces a pair potential 
so that a Q2D mesoscopic superconducting-normal-superconducting (SNS) junction forms 
in the channel. The pair potential is calculated with quasiclassical Green's functions in 
the clean limit. For such a junction alternating Josephson currents and current-voltage 
characteristics (CVCs) are computed, using the non-equilibrium quasiparticle wavefunc- 
tions which solve the time-dependent Bogoliubov-de Gennes Equations. The CVCs ex- 
hibit features found experimentally by the Kroemer group: A steep rise of the current at 
small voltages ( "foot" ) changes at a "corner current" to a much slower increase of current 
with higher voltages, and the zero-bias differential resistance increases with temperature. 
Phase-coherent multiple Andreev reflections and the associated Cooper pair transfers are the 
physical mechanisms responsible for the oscillating Josephson currents and the CVCs. Ad- 
ditional experimental findings not reproduced by the theory require model improvements, 
especially a consideration of the external current leads which should give rise to hybrid 
quasiparticle /collective-mode excitations. 



1. Introduction 

The scattering of electrons into holes and vice versa 
by spatial variations of the superconducting pair poten- 
tial, which occurs at the interfaces between normal (N) 
and superconducting (S) regions Jl], |[ ||, influences sig- 
nificantly the electronic structure and the transport prop- 
erties of mesoscopic superconducting-semiconducting- 
superconducting (SSmS) heterostructures. (In this con- 
text "electrons" and "holes" are quasiparticle excitations 
with momenta above and below the Fermi surface of the 
degenerate semiconductor.) This scattering, also called 
Andreev reflection, establishes phase-coherence between 
the superconducting regions and transfers Cooper pairs 
across the semiconducting layer(s), thus giving rise to 
Josephson currents. 

Most experiments on SSmS junctions involve a quasi- 
twodimensional (Q2D) electron gas in a narrow channel 
between the superconducting banks. The electrons are 
either in the inversion layer which forms between p-InAs 
and Niobium (Nb) (| §, §, or they accumulate in MBE 
grown quantum wells formed by n-type InAs between 
a p-type InAs buffer lay er M , see also ||, or an AlSb 
substrate and two @, |, [Of[jl| or up to 300 0, |l| 
Nb electrodes, separated by insulating AlSb layers. ■ 
Transport properties of SSmS and SNS junctions involv- 
ing high-temperature superconductors are investigated, 
too [§Jj, [d]|. 

One observes peculiar (subgap) structures ||-[@], 
||-[|l2j and spikes in the differential resistances, and 



current-voltage characteristics (CVCs) show first a steep 
rise of the current, then a plateau and finally a linear cur- 
rent increase with "high" voltages |l2|. The observed 
phenomena are interpreted in terms of Andreev reflections 
||-|jjl|. Nitta et al. [fj, ^ assume that the spatially vary- 
ing pair potential, required for these reflections, is induced 
in the InAs by the proximity effect. Chrestin et al. || |(| 
extend the essentially onedimensional, dirty-limit model 
of Aminov et al. to the Q2D case and obtain a gap 
in the density of states which is understood to originate 
from the finite average lifetime of quasiparticles. Volkov 
et al. |l6|] find a small gap in the excitation spectrum of 
the Q2D electron gas, too; multiplying it by the phase dif- 
ference between the superconducting banks they consider 
it as the effective order parameter. 

In this paper we compute the pair potential (order pa- 
rameter) A, induced in the InAs below the superconduct- 
ing Nb electrodes by the proximity effect, from the self- 
consistency equation ([!]). The wavefunctions of quasipar- 
ticles in the Q2D electron gas, Andreev reflected from the 
walls of the superconducting pair potential well, are cal- 
culated from the time-dependent Bogoliubov-de Gennes 
equations. We assume a potential drop eV across the Sm 
region and a model in which the wavefunctions in the Sm 
region evolve from wavefunctions at energies E > A in the 
S banks. With these wavefunctions alternating Josephson 
currents are computed, and their time-averages yield the 
current -voltage characteristics. Comparison with experi- 
ments and a critique of the model terminate the paper. 



2. Proximity Effect 



with the x-dependent critical temperature 



For the numerical computation of the proximity- 
induced pair potential and the current-voltage charac- 
teristics we use system parameters similar to the ones 
reported by the Kroemer group ||, [l(], [l^]. We assume 
that the InAs channel extends in x-direction over 15 nm 
(thickness L x ), and in z-direction the extensions of the in- 
sulating layers between the Nb electrodes, 2a, and of the 
Nb electrodes themselves, D — a, are about 500nm. The 
width L y of the sample in y-direction is of the order of 
magnitude of 100 /im. Furthermore, in the InAs channel 
the "parabolic-equivalent" effective mass of the electrons 
is approximately m* = 0.05 mo (mo: free electron mass) 

We model the electronic structure of these systems in 
the following way. In the Q2D electron gas, below the 
superconducting Nb electrodes, the proximity effect, i.e. 
the diffusion of Cooper pairs, induces a finite supercon- 
ducting pair potential A(z). In the InAs channel, ac- 
cording to Fig. [j]a, the variation of the pair potential 
along the x-axis between the Nb electrode at x = and 
the substrate at x = L x = — 0.08 £o may be neglected. 
(£o = 190 nm is the BCS coherence length of bulk Nb.) 
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Fig. 1 a) The pair potential in a superconducting 
Nb electrode and the proximity-induced pair potential 
in the InAs channel of the samples investigated by the 
Kroemer group |7(j, [7^ /. The critical temperature of Nb 
is Tcs = 9.2 K and that of InAs is Tcn, see text. 
b) The ( Q2D ) SNS-superlattice as a model of the channel 
structure in 

The proximity-induced pair potential in the clean 
limit is calculated with quasiclassical Green's functions 
g(ek,x,oj n ) and the self-consistency equation 



A(x) 



ln(T/T c (x)) + £„l/(n-0.5) 



Tcn 
Tcs 



in N 
in S 



(2) 



The solid angle integration is over all momentum di- 
rections efc. The summation E n over the Matsubara 
frequencies u) n — (2n — Vj'KksT /fi can be cut off at 
hcj n = TiLj c w 10 ksTc (l8|, [n|. n are the Pauli matrices. 
[Eq. ([!]) results from a combination of Kieselmann's |l^] 
eqs. (2.6c) and (2.7b). It has been extended by Bruder 
Jl9[ to the case of anisotropic pairing interactions.] More 



details are given in Ref. 1 20 . We choose the temperature 



T = 2.2 K and use the following material parameters: 
0.21 eV = Fermi energy of InAs, 5.32 eV = Fermi 
energy of Nb, 0.2 = ratio of the Fermi velocity in 
InAs to that in Nb. The temperature-parameter ra- 
tio Tcn/Tcs measures the pairing interaction strength 
in InAs relative to that of Nb. This parame- 
ter is unknown, but we see that even for ratios as 
small as 0.001 one gets a pair potential A in the 
InAs which is about 10% of the Nb pair potential 
Ao = 1.5 meV. Therefore, in the InAs layer supercon- 
ducting and semiconducting regions should alternate in 
z-direction, along which charge transport occurs. As a 
consequence the electronic structure is similar to that of 
ballistic SNS junctions in which only the spatial varia- 
tion of the pair potential breaks translational invariance 
in z-direction, see Fig. lb. 

3. Quasiparticle Dynamics 

The Model 

We use the formalism of the time-dependent 
Bogoliubov-de Gennes Equations (TdBdGE) Q in or- 
der to describe the charge transport in the super-semi- 
superconducting junction formed by the InAs channel 
with the proximity induced pair potential of magnitude 
A. We treat the system as a quasi-twodimensional SNS 
junction. We assume that before times t — t the system 
is in equilibrium, and each quasiparticle state is character- 
ized by a set of quantum numbers k. At times t > i there 
appears a potential difference eV across the semiconduct- 
ing region of length 2a between the superconducting pair 
potentials. It may be due either to a voltage V switched 
on at t = to (voltage bias), or to an applied current which 
increases above the critical Josephson current at t = to 
(current bias). The quasiparticle-wavefunction 



Uk(r,t) 
Vk(f,t) 



(3) 



with electron component (f, t) and hole component 
Vk (f, t) , which evolves from the equilibrium state k un- 
der the influence of the electric field —e z V/2a in the Sm 
region, is determined by the TdBdGE 



(1) 



at 



(4) 



where 



H{r,t) = 



H (r,t) A(z) 
A(z) -H*{r,t) 



(5) 



and 



H (f,t) 



2m* 



7 v + e *4^ 



+ U(x) - pL. (6) 



We model |22|] the spatial variations of the pair potential 
A(z) and of the phase difference <&(t) (which satisfies the 
Josephson equation) by step-functions <d(() (equal to 1, 
if £ > 0, and zero otherwise), using a gauge, where 



A(z) = A-Q(\z\ -a)G(D- \z\) 



is real and 



*(*) 



2eV 



(t - t ) + $ 



Q(a-\z\); 



(7) 



(8) 



$o = ^(*o)j in the case of current bias, is the phase dif- 
ference at which the applied current is equal to the criti- 
cal Josephson current, whereas in the case of a voltage- 
biased junction <& = 0. The scalar potential U(x) de- 
scribes the quantum well which confines the InAs elec- 
trons in x-direction. In an admittedly rather crude 
way we model this confinement by demanding that the 
quasiparticle-wavefunctions vanish in x = and x = L x . 
This means that we assume specular reflection of the 
electrons not only at the band-edge jump (of about 
1.4 eV) from InAs to AlSb but also at the band-edge 
drop (of about 4.9 eV) to the bottom of the Niobium 
conduction band. We assume translational invariance in 
y-direction. The chemical potential /i is that of the en- 
ergy and particle reservoir, to which the channel is cou- 
pled. We neglect all influences of entropy production on 
the chemical potential, because the number of degrees of 
freedom of the reservoir is assumed to be very much larger 
than that of the channel. Then /i is the same as in the 
equilibrium situation and equal to the Fermi energy of the 
InAs channel. 

We expect that the current-voltage characteristics 
calculated for the super-semi-superconducting junction 
modeled by eqs. (f§)-(||) will qualitatively show the same 
features as the ones of superlattices, which may be seen 
as many junctions in series. 

According to the initial condition stated above the so- 
lutions of the TdBdGE (|J) must turn into the station- 
ary quasiparticle-wavefunctions of an SNS junction in 
the limit of vanishing voltage, V — > 0. These wavefunc- 
tions correspond to the ones given in Ref. |2q| . There 
are two classes of stationary states: i) Scattering states 
with absolute value of energy, \E\, larger than the maxi- 
mum value of the pair potential, A, and ii) bound states 
with \E\ < A; E is measured relative to the chemical 
potential fi. The existing theories of ballistic, dissipative 
charge transport in weak links differ in the contribution of 
the two types of voltage-dependent solutions which evolve 
from the two classes of stationary states. In Ref. p4| the 



solutions evolving from the scattering states play the dom- 
inant role, the ones evolving from the bound states matter 
only for very small voltages. (Refs. ^7) and re- 

lated theories consider only the evolution from scattering 
states.) In Ref. |%| , on the other hand, current -voltage 
characteristics (with negative differential resistance ob- 
served in short metallic weak links used for millimeter- 
wave mixers |2^]) were calculated from solutions which 
evolve from the bound states. The physical difference 
between the two classes of theories consists in the as- 
sumption about the rate at which quasiparticles, in a 
sequence of relaxation cycles, start their motion in the 
electric field from energies \E\ < A: In Ref. j24j] quasi- 
particles start after each relaxation cycle, whereas in Ref. 
p2[ they start after each Andreev reflection. The idea 
behind the assumption of Ref. |^| is the following: If the 
SNS junction is connected to the current or voltage source 
by normal conducting external current leads, practically 
all quasiparticle excitations, forming the shifted Fermi 
sphere of the current carrying normal leads, have ener- 
gies \E\ < A so that they are totally Andreev reflected at 
the external interfaces between leads and junction. When 
decaying in the superconducting banks of the junction 
they form Cooper pairs and induce collective (supercur- 
rent carrying) modes [STL [30|. These collective modes 
should turn again into quasiparticle excitations within the 
N layer of the junction, with initial energies \E\ < A as 
in the external leads. Because of charge conservation the 
rate at which the collective modes turn into quasiparticle 
excitations is equal to the rate of Andreev scattering at 
any of the interfaces. 

The BCS theory of superconductivity and the related 
TdBdGE describe a superconductor as a grand canonical 
ensemble. If the Cooper pairs created (destroyed) by An- 
dreev reflections in an SNS junction directly go to (come 
from) a reservoir of Cooper pairs, with no normal current 
leads between junction and reservoir, the model of Ref. 
PH should apply. If, on the other hand, the reservoir is 
connected to the junction by normal conducting current 
leads, the model of Ref. |22j should be appropriate. In the 
latter case one should analyze in detail the spatial distri- 
bution of the hybrid excitations, consisting of quasiparti- 
cles (in the normal current leads and the normal region of 
the junction) and collective supercurrent carrying modes 
(in the S layers). This is being investigated presently. In 
this paper we disregard, like Refs. ^4|-|27|], the possible 
influence of normal current leads and assume that only 
quasiparticles with energies \E\ > A are incident from 
the superconducting banks onto the NS interfaces of the 
junction. 



Multiple Andreev Reflections 

There are four different types of incident quasiparticles 
with energy E > A. Their wavefunctions, which satisfy 
eq. (^) in the superconducting banks and which are nor- 
malized to the volume L x L y (D — a) of one S layer, are 



given by The solutions of eq. (J4|) in the normal layer with the con- 

stant electric field are 



r, t) = V (x, y) ■ e'^N^E) 
x e m[k *f + P^ E)]z O(-af3z - a)Q(D - \z\) (9) 



*f(E;r,t) = r,(x,y) ■ e -^ E+ ^ v ^C^ {E) C^, 



with 

77(1,1/) = sin(k x x)e' tk « y ; k x = — s, s = 1, 2 , (10) xe L ^ J e' 1 ^ *° « 6(a - \z\) . (17) 

describing the plane waves in y-direction and the stand- They result from an expansion of Airy functions ||l| in 
ing waves in the quantum well forming the channel. At the limit 
low temperatures, for /j — 0.2 eV, L x = 15 nm, and h 2 k 2 

2m 



to* = 0.053 too, only the subband states with s = 1 and IT^T > \E±eV\, A . (18) 



s = 2 are occupied. The indices a and /3 in eq. (|9j) indi- 
cate the character of the quasiparticle: it is electron-like The coefficients Da p {E) and {E) are determined 

[P = +1) or hole-like {fi = -1) with positive (a = +1) by matching the superp osition 
or negative (a = —1) momentum in z-direction. Thus, 

the complete set of quantum numbers characterizing the oo 

quasiparticle-wavefunction (|) is k = {s;k y ; {E,a,/3)}. $g(r,t)=^ / dE j*^ {E:r,t) + ^f t (E;r,t)\ 

1 p-L 

N "W = \ T-rm — utt rwrn -Q(^- a ) a 1 ) « ^ ( 19 ) 

y L x Ly (D - a) [1 + ^ formed from the solutions (§) and (|16|), to the superposi- 



is the normalization factor, with 



tion 



7(E) _ b-w^j [or |E| , A . (12) nc.o-x; / 



(20) 



For < A formed from the solutions (0), at the NS interfaces at 

, z = ±a: 

7(S) - • (13) 

*g(* ? ,*)i,=±. = n(*'.*)u ± a ■ ( 2l ) 

7(-E) is the probability amplitude of Andreev reflection 

of a quasiparticle with energy E. [Arguments of ~f(E) These matching conditions correspond to the ones used 
less than A appear in the probability amplitudes of in Refs. |2^, |3lj]. 

multiple Andreev reflections A^(E), see eq. (p7|).] The The integrals over E also contain the solutions which 

wave number of quasiparticle propagation in z-direction evolve from negative energies — \E\. The corresponding 
a [k z f + /3k(E)] is determined by the component of the wavefunctions 
Fermi momentum perpendicular to the interfaces 

a fi , IP1 .,._K-.(-| £ li r >*) 



k zf = JkJ - fc 2 (14) *S%HE\;r,t) = 



and 



are obtained from the ones of positive energy +|2£|, 



s i gn (£)V|!EAl for |£|> A 
K{E) = \ .;43! "" taH<A ■ <15) ^S5.<^'..) 



— ir~ 
m 



where = fi, k e = e x k x + e y k y , and v z f = -^j-. 

The wavefunctions for the quasiparticles which are An- by the relation 
dreev reflected at the NS interfaces and travel back in the 

S banks are: (u a s %{-\E\;r, t)\ ( v~ s ^ *(+\E\: 7. 1 , 

r, t) = V(x, y) ■ e-* Et D"P(E) 



1 \ \v^ m (-\E\;r,t)J \-u s % f3 *{+\E\^t)) ' ^ 

The negative energy states are the ones occupied in the 
x e la[k 'f +ME)]z Q(aPz - a)G(D - \z\) . (16) ground state of the system. 



With the wavefunctions (|9|) , (|16| ) , and (|17j) the match- 
ing conditions (^l|) become: 

oo 

a J \i\ ) J 



x {N {E)Q(-a(3z-a) + D af3 (E)Q{af3z-a)}\ z=±a 



a/ J V 9 



x e 



*.H h-r~. — 



e -4/3'*of e(a _ | z | 



(23) 



We do not need to know the coefficients D a @(E) explicitly, 
because we calculate the current density in the normal 
layer of the junction. After straightforward but tedious 
algebra one finds from eq. (p3|): 



<?# = t£#? ( E -[2n + l- \m°*V 



n=0 



(24) 



with 



cfe'(E)=cfP'(E)e 



x e 



A" E 



(25) 



Here 



zf P '(E) = N\E) [1 - j(E) 2 ] e -i*W* e -ip-mi°*o 



x e 



(2- J 3/3')(E-ia/3/3'eV) + |c,eV 



7(S + aey) 



(26) 



and 



l;;(£) ^ n7(^ + K|K) fornX) ^ 
1 for ri = 

is the multiple Andreev reflection probability amplitude 

(7^' (i?) consists of the phase-coherent superpositions 
of all amplitudes df?' (E - [2n + 1 - \l3f3']aeV) which 
originate from incident quasiparticles described by the 
wavefunctions - [2n + 1 - \j3(i']aeV; f, t) . In 

other words: The probability amplitude C a @ (E) of find- 
ing a quasiparticle (characterized by a and /?') in the 
N layer, with energy in the range between 
E — eV/2 and E + eV/2, is the sum of the ampli- 
tudes c^ 1313 ' (E - [2n + 1 - \(3f3')aeV) ; |c«^'| 2 gives 



the probability that a quasiparticle, incident from one 
of the S-banks and characterized by a, (3, and energy 
E — [2n + 1 — ^P(3']aeV, reaches the energy range be- 
tween E — eV/2 and E + eV/2 in the N layer after n 
Andreev-reflection cycles, being an electron (hole) for 
(3> = +l(-l). 

However, the most convenient way of calculating 
the current density f in the N layer is not the one 
which uses the wavefunctions with C a/3 (E). It is much 
easier to count properly all quasiparticle contributions 
to j of cq. (||^) if one uses new N layer wavefunc- 
tions W Nk (r, t) . (These are the wavefunctions (^) in the 
N layer.) They are obtained by adding up the phase- 
coherent wavefunction components which, by multi- 
ple Andreev reflections, originate in the N layer from 
one particular incident quasiparticle described by eq. @ 
with a, (3, and energy Eq- The sum of these components 
is given by the right-hand side of eq. (|l7]) if there one sets 
E = E + [2n + 1 - \(3f3']aeV and replaces C a/3 ' (E) by 
the right-hand side of eq. (24) with one change: Instead 



of summing over the index fi characterizing the nature 
(electron-like or hole-like) of the incident quasiparticle 
one has to sum over the index (5' which characterizes the 
quasiparticle nature (electron or hole) in the N layer. 
One can show: i) The new wavefunctions 

oo 

9 Nk (f>, t) = V (x, V )-J2 ^ e -i[ £ o + (2n+l-i/3/3')^+i/3'e^] t 



n=0 13' 



- 1±£ \ la 



xc^ (E ) ,4 e 



1 eV j.+/3'[s + (2»+l- J/3/3')<*ev] 



x e-^ /3 '*°»e(a- \i 



(28) 



match to the solutions in the S layers, ii) In the limit of 
vanishing voltage, V — > 0, the absolute squares of the elec- 
tron and hole components of ^Nk exhibit the resonances 
of the virtually bound states one has in the stationary 
situation 

4. Current Voltage Characteristics 

In the Bogoliubov-de Gennes-formalism the current 
density in the N layer is : 

3 t) = — ^2 {fk u Nk(r, t)p e u N k(r, t) 

m k 



[1 - fk] v N k (f, t)p e v% k (r , t) + c.c. } 



(29) 



where the sum over k = {s; k y ; ( E=Eq, a, /?)}, character- 
izing the incident quasiparticles, only counts the states 
with positive energy E > A and one spin direction; 

fk = ^exp(-^^) + lj . Applying the momentum op- 
erator p e = 5V + e z ^$(t) to the wavefunctions given by 



cq. fl28|), averaging over the channel thickness L x , and 
defining 



1 



VE 2 - A 2 l + 7(^) 2 



7(£ + aeV) 



x Re ^ ( 

n,n'— 



.0/ 2(n-n') e Vm* . 



-m(n — n )<&o 



Xc 



4(n-n') E+ („ + „' + l- l^^'jceV 



x^,*(£- -a/3/3') ■ a 



A«(£--a/3/3'ey) 



%xlk 2 F -kl 



+ 



_h_ lh.2 _ 1,2 



@(a-\z\) 



we obtain the current density in the N layer as 



(30) 



T(z,t) 



E 

a/3 



dE 



2 

E 

s=l 



dA y (fJ^XE-^t) 



(31) 



As usual we integrate over all energies up to the Debye 
energy huj D ; 



h = \ h 2 — 2™ 

""y,max — \ n p ^: 



The time average of j(z, t) is 



■* M 

= ^ J dtj(z,t) 



(32) 
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Fig. 2 The time- and voltage-dependent current density 
at temperature T = 2.2 K in the quasi-twodimensional 
electron gas of an SSmS junction with the proximity- 
induced pair potential of Fig. The length of the 
Sm-layer is 2a — 500 nm, and the electron density is 
5-10 ls cm- 3 . 
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Fig. 3 The time- and voltage-dependent current density 
at temperature T — 2.2 K in a threedimensional SNS junc- 
tion. The length of the N layer is 2a = 500 nm, and the 
Fermi energy is 0.21 eV. 



with Tm —>■ oo. In numerical calculations we choose 
Tm = 1 s, which is very much larger than all charac- 
teristic time scales. 

The time-dependent current density J{0, t) in the 
center of the junction is computed by evaluating nu- 
merically equation d|l]) for the voltage-bias case, i.e. 
$0 = 0- We assume a proximity-induced pair potential 
A = 0.3 meV. The other parameters are indicated in 
Section 2. The result is shown in Fig. |[ For com- 
parison Fig. H shows the time-dependent current den- 
sity of a three-dimensional SNS junction with the same 
pair potential but translational invariance in x— and y— 
directions |3^] . Its basis is the 3D current density equa- 
tion which one obtains from eq. ( J3l| ) essentially by 
changing the sum over the indices s=l and 2 of the 
two occupied subbands into an integral over continous 
fc^-values, and by replacing the effective mass m* by the 



free electron mass mo- Not surprisingly, there are only 
quantitative differences: The fundamental mechanism of 
multiple Andreev reflections, responsible for the charge 
transport in mesoscopic weak links, works the same way 
in three- and quasi-twodimensional electron gases. The 
amplitudes of the current density oscillations over time 
are small. Therefore, the time-averaged current densities 
j{0), computed from eq. ( |32|) and presented in Fig. ||, 
result in current -voltage characteristics (CVCs) which do 
not differ much from the CVCs in Fig. || at any fixed 
time. (In order to save computer time, the CVCs of Figs. 
H and H have been calculated with less numerical accuracy 
than those of Fig. [|. Therefore, the latter are smoother 
than the former. Also, the low-voltage slope in Fig. || is 
not as steep as in Fig. || because only four, instead of 36, 
Andreev reflection cycles have been taken into account in 
its computation.) 
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Fig. 4 Current-voltage 
SSmS junction at T 



of 



Q2D 



0.4 0.5 0.6 
V[mV] 

characteristics 
= 2.2 K (lower curve) and 
T = K (upper curve), same parameters as in 
Fig. [|. The inset shows the steep rise of the current 
density with small voltages; again, the lower curve is for 
the higher temperature. 

5. Discussion 

The computed current densities in Figs. || - | agree 
with the experimental findings of Thomas et al. [[L2| in 
the following qualitative aspects: i) There are voltage- 
dependent current-oscillations in time. The period is the 
inverse of the Josephson frequency v = 2eV/h: The dom- 
inant contributions of the multiple Andreev reflections to 
the alternating (Josephson) current come from the terms 
with n — n' = 1 in eq. (pQ). ii) The CVCs of the time- 
averaged current densities exhibit what Thomas et al. call 
the "corner current" . This current value separates the 
steep, nearly linear increase of the current with very small 
voltages (also called the "foot" in prior investigations of 
non-equilibrium effects in weak links |22 ) from a wide 
range of slow current increase with higher voltages. The 
physical mechanism responsible for the drastic change of 
current growth with voltage is the decreasing number of 
relevant Andreev reflection cycles in the energy range be- 
tween —A and +A with increasing voltages above the 
voltage of the "corner current" . iii) The zero-bias dif- 
ferential resistance dV/dI\v=o increases (and the current 
decreases) with temperature. These points of qualitative 
agreement between experiment and theory support the 
view of Thomas et al. that multiple Andreev reflections 
are responsible for the observed noncquilibrium effects 
(12). The characteristic shape of the CVC described in 
ii) has been also observed for a single SNS junction by 
Kroemer et al. fic| . There is even satisfactory quanti- 
tative agreement between the magnitude of the experi- 
mental current densities and voltages in Fig. 2 of Ref. 
|To| and the computed ones in our Fig. |J. Not surpris- 
ingly, quantitative agreement with the measurements of 
Thomas et al. performed on more than 300 SNS junctions 
in series Jl2| is not so good: Our current densities are at 
least four times greater than theirs. For many SNS junc- 



tions in series the following approximations of our model 
may have to be modified: It a) confines all electrons to 
the channel, neglecting their penetration into the Nb con- 
tacts, b) neglects trapping of electrons by surface states 
in y = and y — L y , and c) disregards the small proba- 
bility of conventional scattering from spatial variations of 
the pair potential |33|, |^] , which may occasionally disrupt 
the Andreev reflection cycles |T^] . 

Two important qualitative features of the experimen- 
tal CVCs are not reproduced by our model. The first is 
the dependence of the steep slope of the CVCs for van- 
ishing voltages (the "apparent activation energy" in the 
Arrhenius-plot language of Thomas et al.) on the sam- 
ple width L y . In order to include effects of this width 
in the model, one may have to take into account - via a 
phenomenological, Lj,-dependent relaxation time - that 
the phase-coherence of Andreev reflections can be de- 
stroyed by the temporary trapping of electrons in sur- 
face states. Alternatively, one may also think of inhomo- 
geneities in the current density, caused by the (neglected) 
magnetic field of the current and the finite Josephson pen- 
etration depth. The required improvements would not 
change the fundamental character of the model. A signif- 
icant model change, however, may be necessary in order 
to reproduce the second characteristic feature of the ex- 
perimental CVCs not present in our figures: The nearly 
flat plateau in the CVCs between the "corner current" at 
the end of the low-voltage "foot" and the linearly rising 
excess current at high voltages. Such a plateau is ab- 
sent in all theories where the principal contributions to 
the Multiple- Andreev-Reflection current are from quasi- 
particles which evolve from (scattering) states with en- 
ergies above the gap. Our CVCs, calculated with solu- 
tions of the time-dependent Bogoliubov-de Gennes Equa- 
tion, and that of Gunsenheimer and Zaikin, obtained 
with the Keldysh technique [Q, are qualitatively quite 
similar and lack the plateau. On the other hand, CVCs 
computed with quasiparticle wavefunctions which, at the 
rate of Andreev reflection occurrence, evolve from states 
with energies below A show the plateau for appropriate 
inelastic mean free paths ^2|. As discussed in Section 3, 
this should correspond to a situation where the sample 
is connected to the external energy and particle reser- 
voir by normal conducting current leads, and where hy- 
brid quasiparticle/collective-mode excitations at energies 
smaller than A exist. The inclusion of the subgap start- 
ing states in the theory, also necessary in order to obtain 
the direct Josephson current in the case of current bias, 
$0 7^ 0, and vanishing voltage, is a subject of work in 
progress. 
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